Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation          下载免费PDF全文

Tongbi Tu  Ali Ercan  M. Levent Kavvas 《水文研究》2018,32(10):1406-1419

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.  相似文献   

The Upper Devonian Kellwasser event recorded in a regressive sequence from inner shelf to lagoonal pond,Catalan Coastal Ranges,Spain     

《Sedimentology》2018,65(6):2055-2087

This study, conducted in the Catalan Coastal Ranges, north‐east Spain, describes the Upper Devonian Kellwasser event in a shallowing‐upward sequence of black shales, siltstones and quartz arenites. This sequence was deposited in a progradational and regressive coastal system where the sedimentary environment evolved from the inner shelf to a lagoonal pond located landward of the shoreline. Three anomalous succeeding steps have been identified by geochemical analysis. The first one, detected on the inner shelf, was characterized by oxygen depletion and high organic productivity. The second, detected in the nearshore, was caused by hydrothermal activity occurring under normal oxic conditions. The third and most intense step was identified in the muds of the lagoonal pond and has been linked to strong anoxic conditions, elevated clastic input derived from changes in the weathering regime at the source area and moderate hydrothermal activity. The Kellwasser event is thus defined in the study area as stepwise and multi‐causal. This is the first time that the Kellwasser event has been identified in a sedimentary environment behind the shoreline. It is also the first time that it has been reported in the Catalan Coastal Ranges.  相似文献   

Three dimensional simulation of fluid flow in X‐ray CT images of porous media     

A. Al‐Omari  E. Masad 《国际地质力学数值与分析法杂志》2004,28(13):1327-1360

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction     

M. Zeghal  U. El Shamy 《国际地质力学数值与分析法杂志》2004,28(14):1361-1383

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

Shearing flows of a dry granular material—hypoplastic constitutive theory and numerical simulations     

Chung Fang  Yongqi Wang  Kolumban Hutter 《国际地质力学数值与分析法杂志》2006,30(14):1409-1437

  相似文献   

Luminescence dating of sand–loess sequences and response of Mu Us and Otindag sand fields (north China) to climatic changes     

Yali Zhou  Huayu Lu  Jiafu Zhang  Joseph A. Mason  Liping Zhou 《第四纪科学杂志》2009,24(4):336-344

The sand–loess transition zone in north China is sensitive to climate change, and is an ideal place to investigate past environmental changes. However, past climate change at millennial–centennial timescales in this region has not been well reconstructed because of limited numerical dating. Alternations of sandy loam soils with aeolian sand layers in the Mu Us and Otindag sand fields, which lie along the sand–loess transition zone, indicate multiple intervals of dune activity and stability. This change is probably a response to variations of the East Asian monsoon climate during the late Quaternary. The single aliquot regeneration (SAR) optically stimulated luminescence (OSL) dating protocol, which has been successfully applied to aeolian deposits worldwide, is applied to these two sand fields in this study. The OSL ages provide reliable constraints for reconstruction of past climate changes at suborbital timescale. Sections in both sand fields contain aeolian sand beds recording millennial‐scale episodes of dry climate and widespread dune activation, including episodes at about the same time as Heinrich Event 5 and the Younger Dryas in the North Atlantic region. These results demonstrate the potential of aeolian sediments in semi‐arid north China to record millennial‐scale climatic events, and also suggest that dry–wet climate variation at the desert margin in China may be linked to climatic change elsewhere in the Northern Hemisphere, through atmospheric circulation. This article was published online on 27 November 2008. An error was subsequently identified. This notice is included in the online and print versions to indicate that both have been corrected (16 December 2008). Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

Biostratigraphic and aminostratigraphic constraints on the age of the Middle Pleistocene glacial succession in north Norfolk,UK     

Richard C. Preece  Simon A. Parfitt  G. Russell Coope  Kirsty E. H. Penkman  Philippe Ponel  John E. Whittaker 《第四纪科学杂志》2009,24(6):557-580

Considerable debate surrounds the age of the Middle Pleistocene glacial succession in East Anglia following some recent stratigraphical reinterpretations. Resolution of the stratigraphy here is important since it not only concerns the glacial history of the region but also has a bearing on our understanding of the earliest human occupation of north‐western Europe. The orthodox consensus that all the tills were emplaced during the Anglian (Marine Isotope Stage (MIS) 12) has recently been challenged by a view assigning each major till to a different glacial stage, before, during and after MIS 12. Between Trimingham and Sidestrand on the north Norfolk coast, datable organic sediments occur immediately below and above the glacial succession. The oldest glacial deposit (Happisburgh Till) directly overlies the ‘Sidestrand Unio‐bed’, here defined as the Sidestrand Hall Member of the Cromer Forest‐bed Formation. Dating of these sediments therefore has a bearing on the maximum age of the glacial sequence. This paper reviews the palaeobotany and describes the faunal assemblages recovered from the Sidestrand Unio‐bed, which accumulated in a fluvial environment in a fully temperate climate with regional deciduous woodland. There are indications from the ostracods for weakly brackish conditions. Significant differences are apparent between the Sidestrand assemblages and those from West Runton, the type site of the Cromerian Stage. These differences do not result from contrasting facies or taphonomy but reflect warmer palaeotemperatures at Sidestrand and a much younger age. This conclusion is suggested by the higher proportion of thermophiles at Sidestrand and the occurrence of a water vole with unrooted molars (Arvicola) rather than its ancestor Mimomys savini with rooted molars. Amino acid racemisation data also indicate that Sidestrand is significantly younger than West Runton. These data further highlight the stratigraphical complexity of the ‘Cromerian Complex’ and support the conventional view that the Happisburgh Till was emplaced during the Anglian rather than the recently advanced view that it dates from MIS 16. Moreover, new evidence from the Trimingham lake bed (Sidestrand Cliff Formation) above the youngest glacial outwash sediments (Briton's Lane Formation) indicates that they also accumulated during a Middle Pleistocene interglacial – probably MIS 11. All of this evidence is consistent with a short chronology placing the glacial deposits within MIS 12, rather than invoking multiple episodes of glaciation envisaged in the ‘new glacial stratigraphy’ during MIS 16, 12, 10 and 6. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

A global wave parameter database for geophysical applications. Part 1: Wave-current–turbulence interaction parameters for the open ocean based on traditional parameterizations     

Nicolas Rascle  Fabrice Ardhuin  Pierre Queffeulou  Denis Croiz-Fillon 《Ocean Modelling》2008,25(3-4):154-171

Ocean surface mixing and drift are influenced by the mixed layer depth, buoyancy fluxes and currents below the mixed layer. Drift and mixing are also functions of the surface Stokes drift U_ss, volume Stokes transport T_S, a wave breaking height scale H_swg, and the flux of energy from waves to ocean turbulence Φ_oc. Here we describe a global database of these parameters, estimated from a well-validated numerical wave model, that uses traditional forms of the wave generation and dissipation parameterizations, and covers the years 2003–2007. Compared to previous studies, the present work has the advantage of being consistent with the known physical processes that regulate the wave field and the air–sea fluxes, and also consistent with a very large number of in situ and satellite observations of wave parameters. Consequently, some of our estimates differ significantly from previous estimates. In particular, we find that the mean global integral of Φ_oc is 68 TW, and the yearly mean value of T_S is typically 10–30% of the Ekman transport, except in well-defined regions where it can reach 60%. We also have refined our previous estimates of U_ss by using a better treatment of the high frequency part of the wave spectrum. In the open ocean, U_ss  0.013U₁₀, where U₁₀ is the wind speed at 10 m height.  相似文献   

The Upper Pleistocene to Holocene sediments on the Mediterranean island of Lampedusa (Italy)     

C. Giraudi 《第四纪科学杂志》2004,19(6):537-545

The island of Lampedusa lies on the northern edge of the African continental shelf, but during some Quaternary marine lowstands it was joined to the African continent. The study and dating of the aeolian, alluvial, detrital sediments, calcareous crusts and speleothems have established that the climatic–environmental variations recorded on the island can be related chronologically to those known for northern Libya, Tunisia and the Italian peninsula. During the Last Glacial Maximum, phases of Saharan dust accumulation on Lampedusa occurred, and were coeval with dust accumulation in crater lakes and on high mountains in central‐southern Italy, and with phases of glacial advance in the Apennines and in the Alps. During the late Holocene, accumulation of Saharan dust on Lampedusa occurred but there was little accumulation of dust on the northern side of the Mediterranean Sea. With the new data from Lampedusa, it is possible to envisage two different scenarios of atmospheric circulation relating to the Last Glacial Maximum and to the late Holocene. During the Last Glacial Maximum, southerly atmospheric circulation brought rainfall to the southern slopes of the Alps and to the Apennines. During the late Holocene, a prevalent westerly atmospheric circulation became established in the northern Mediterranean. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

The Early Pleistocene Praetiglian and Ludhamian pollen stages in the North Sea Basin and their relationship to the marine isotope record     

T. Meijer  P. Cleveringa  D. K. Munsterman  R. M. C. H. Verreussel 《第四纪科学杂志》2006,21(3):307-310

In this paper, the currently accepted correlation of the Early Pleistocene Ludhamian stage of England with the Tiglian‐A sub‐stage of the Netherlands is challenged. Recent investigations of Early Pleistocene marine North Sea deposits from a borehole near Noordwijk (the Netherlands) yielded evidence from molluscs, dinoflagellate cysts and sporomorphs for an alternation of warm‐temperate and arctic intervals within the Praetiglian and Tiglian stages. Marine equivalents of the terrestrial‐based pollen sub‐stages Tiglian A and B have been recognised in the upper part of the sequence. A Praetiglian age can be assigned to the lower part of the sequence on the basis of mollusc analysis. Within the Praetiglian, an alternation of warm and cold phases has been recognised from both the dinoflagellate cyst and molluscan records. Three cold phases within the Praetiglian are tentatively correlated with marine isotope stages (MIS) 96–100. The molluscan assemblages provide evidence for climate forcing of the sea level: highest sea levels are reached in the warm‐temperate intervals. Within the Praetiglian, an interval with an acme zone of the dinoflagellate cyst Impagidinium multiplexum, is correlated with the Ludhamian and tentatively linked to MIS 97 and/or MIS 96. The cold molluscan assemblages from the Noordwijk borehole include an acme zone of Megayoldia thraciaeformis, the first and only occurrence of this North Pacific bivalve in the North Sea Basin. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献